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We calculate the cross sections for reactions of the J/ip with light mesons. We also evaluate 
its finite temperature spectral function. We investigate separately the role of elastic and inelastic 
channels and we compare their respective importance. We describe J/ip absorption channels that 
have not been considered previously to our knowledge. The relevance of our study to heavy ion 
■ collisions is discussed. 
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The study of relativistic heavy ion collisions offers the tantalizing possibility of observing many-body effects in a 
strongly interacting system at densities and temperatures far removed from equilibrium. At ultrarelativistic energies, 
the main focus of the active experimental and theoretical programs is the creation, observation, and interpretation of 
. a new form of matter: the quark-gluon plasma. Its existence is a prediction of QCD, even though some ambiguities 
concerning the specific nature of an eventual phase transition and its experimental signatures Q still remain. It is 
fair to say that the activity generated by this field makes it one of the most exciting areas of contemporary subatomic 
physics. 

The J/ip meson has been singled-out as a promising candidate to signal deconfinement. Indeed, the presence of a 
high temperature quark-gluon plasma would screen the cc interaction B or ionize the quarkonium state J3| , leading 
to a suppression of the J/ip in events where the plasma is produced in comparison with events where it is not. While 
the suppression of J/ip (and ip') in p-A and in heavy ion collisions involving medium-mass projectiles at 200A GeV 
can be explained by absorption models without plasma assumptions (|||, the subsequent Pb + Pb data has led 
to analyses involving plasma formation However, the plasma interpretation of this "anomalous" J/ip absorption 
observed with the Pb projectile needs the introduction of model-dependent assumptions. Furthermore, alternative 
explanations which rest purely on hadronic grounds are starting to appear Any model of J/ip suppression, 

whether it includes hadronic "comovers" || or not, relies on simple assumptions about the size of the cross sections 
with nucleons and light hadrons. Unfortunately, up to recently only a few calculations of the interaction cross section 
of J/ip with light hadrons could be found in the published literature [ ^o|JTT| , and the results of those calculations are 
not in agreement with each other ]T^ |. 

Our aim in this work is the following. We plan to systematically explore the different channels of interaction of the 
1 J/ip with light hadrons and calculate the corresponding cross sections, in light of recent calculations with an effective 
hadronic Lagrangian (l^Jl^]. We shall perform no attempts to find heavy ion data here. However we will try to bring 
the study of the hadronic interactions of the J/ip meson closer to the level of sophistication that the lighter vector 
mesons (p, u, and <p) currently enjoy. Consequently, we calculate the spectral function for the charmonium state, in 
a gas of light mesons at finite temperature. Bear in mind that this does not imply that the J/ip is thermalized. We 
view this calculation as a necessary prelude to a more complete understanding of the behavior of the charmonium 
bound states in hadronic matter at finite temperature and density. Because of the direct decay into muon pairs, the 
J/ip spectral function is directly accessible experimentally. 

Our paper is organized as follows: in the next section we discuss details of a heavy meson chiral Lagrangian bearing 
hadronic interactions upon which our quantitative estimates are based. We also describe a slight variant of this model 
that has been used in phenomenological applications. The dominant channels in this study, both elastic and inelastic, 
are then considered and the associated cross sections are shown. We introduce new channels, to our knowledge, for 
J/ip absorption on mesons. We then proceed to a discussion of the scattering widths induced by the interactions. 
We will then show the resulting J/ip spectral function, and explore its temperature and momentum dependence. We 
summarize and conclude. 
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II. CHIRAL LAGRANGIAN FOR LIGHT AND HEAVY MESONS 



We discuss here the basic assumptions and ingredients in our chiral Lagrangian approach for light and heavy 
pseudoscalar and vector mesons. We shall model the interaction of the charmonium state with lighter mesons through 
meson exchanges. In order to include charmed mesons, the smallest possible symmetry group that potentially contains 
the relevant phenomenology is SU(4). However, SU(4) is in fact badly broken by the large mass of the charmed quark. 
This also can be seen in the poor agreement obtained between the extended mass formula and the experimentally 
measured masses [ fL5[ . We adopt here the following pragmatic viewpoint: we work here with the physical mass 
eigenstates and the physical mass matrix will represent the relevant breaking of the original symmetry. Furthermore, 
we compare two calibration methods, (i) the chiral gauge coupling will be uniquely determined from light vector 
spectroscopy, namely the p meson, and (ii) relevant coupling constants are individually chosen by either empirical 
constraints where they exist or model calculations in the absence of measurement. We feel that it is crucial for this 
effective approach to be in tune with the largest possible range of phenomenology at the appropriate energy scale. 

Description of light and heavy pseudoscalars in a single framework can be obtained using a four-flavor chiral 
Lagrangian. The basic nonlinear SU(4) er model apart from mass terms is 



£o = -^■Tr(d li U& t U*), 
U = exp — 



(2.1) 



The constant F n ~ 135 MeV and (j> is the pseudoscalar multiplet matrix. Correct normalization leads to 
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d„U -> V^U = d„U - igA L U + igUA R , 



and we add kinetic terms 
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where and A R are the chiral spin-1 fields and where 



F RpR^ +7Tr (F^UF^tf) , 



(2.3) 



(2.4) 



F^ = d lt Ai-a v A^-ig[A^,Ai] , 
F R =d^A R -d u A R -ig [A*, A*] 



Next we add mass terms for the spin-1 fields plus two generalized mass terms 
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The aim for the present model is to describe the normal parity states, so we must eliminate the axial- vector matrix 
field A^ = A^ — A R . To accomplish this, we follow the ideas presented in Ref. |ll| and make a gauge transformation 
resulting in A' — 0. Equivalently, in the primed gauge, A^' = A Rl = p^. The vector meson matrix multiplet is 
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The specific choices C/ 1 / 2 = £, f7 1 / 2 = £t and 
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4 = CM 1 + 

a* = ^ + 

U = S1£ (2.8) 



will "gauge away" the positive-parity states from the model by producing the requisite A'^ — 0. 

Gauging away the axial fields with the above mentioned transformation yields the following Lagrangian (utilizing 
Hermiticity <f> = (fy and p p = p' ) 

C Q = 

C 1 = ( 7 -l)Tr[F fll/ (p)F^(p)}, 

£ 2 = (B + 2C- 2ml) Tr (p^) + 2^-20-2™*) ^ ^ ^ ^ 

+ ^Tr([^ p<f ) - {B + 2 ° + 2m2o) Tr (0, W ) , (2.9) 



where we have defined 



F»u(p) = d p p„ - d v p p - ig [p p , p v \ . (2.10) 



We notice that under the gauge transformation the original kinetic piece for the pseudoscalars vanishes — it reappears 
in £ 2 - >Ci becomes a kinetic Yang-Mills term for the field p p , while its correct normalization points to 7 = 3/4. 
Furthermore, £2 includes the mass terms for the p p field, kinetic terms for 0, as well, it includes three-point and 
four-point interaction terms. It is clear that the mass term for p p spoils local gauge invariance. In order to leave 
Co + Ci + £2 chirally and locally gauge invariant, we choose B + 2C — 2ml = 0- Omitting kinetic energy terms, we 
arrive at the model's chiral and gauge invariant set of interactions. They are the following 

Ant = igTcfa r<M) - ^Tr([0,p^] 2 ) +igTt(a„p v [p^p v ]) + ^-Tr{[(f , p v f) (2.11) 

where, for convenience, we have attached physical significance to the gauge coupling constant through 

2{B-2C-2ml) _ g /mrr 
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= 9- (2.12) 



Before leaving the formalism we summarize up to this point. An effective chiral Lagrangian of light plus heavy 
pseudoscalar and vector mesons has been constructed to be fully chirally U(4)xU(4) invariant. In particular, no loose 
ends are left in the model since the axial fields have been completely gauged away. The result is w ritten in Eq. ([2.9]) . 



From here we further impose local gauge invariance and arrive at interaction terms given in Eq. (2.11). The model 
has two input parameters: F T and g p7r7r . Since we adjust (7 p7r7r using the decay rate into pion pairs at the physical rho 
pole, we implicitly also use m p . The matrix algebra implied in the above expressions can now be explicitly carried 
out to obtain Lagrangians involving specific physical fields. 



A. Chiral Model Predictions 



The model is constructed in some sense to pivot off the rho meson since p decay into two pions uniquely determines 
the gauge coupling constant g. Ta king T p — 151 MeV, and m p — 770 MeV, we find g P w — 2g ~ 8.54. Having fixed 
the single parameter in Eq. ( 2.11 ), we are in position to "predict" widths for meson decays where phase space is 
open, and we are particularly interested in the strangeness and charm sectors. K*'s and D*'s are chosen to test the 
symmetry breaking effects in the extremes. Results for widths are listed in Table |. We find K* widths consistent 
to within 10% of experiment, and D* widths consistent with other model calculations [|i7|-|l9|. We are hesitant to 
read too much into these numbers, but they begin to suggest that the symmetry breaking effects might largely be 
accounted for merely by using the physical mass eigenstates. 

All coupling strengths are now in principle fixed. We have first evaluated absorption cross sections for reactions 
involving the J/tp in the initial state. We list the processes under consideration in Table ||. Note that each process 
can in principle involve several Feynman diagrams and their interference. 
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For the sake of brevity, we do not show all the cross sections we have evaluated. We shall restrict ourselves 
here to the most important ones. Our findings first support the notion that the elastic channels are quantitatively 
small. Apart from the ui to be discussed later, the largest contribution in this category is shown in Fig. |], and is 
J/ip + p — > J/ip + p. By y/s = 6 GeV, the cross section has just about risen to f mb. This process is modeled 
by pion exchange and therefore involves vector-vector-pseudoscalar interactions which are not included in the chiral 
Lagrangian. Instead, the relevant Lagrangian is 



C = ge a ^ v d a V p d^V lJ <P . 
The coupling strength is fixed to match the measured (26) J/ip — > irp decay width. 



(2.13) 
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FIG. 1. Elastic cross section J/ip + p — ► J/ip + p. 



The elastic process immediately following this one in importance is J/ip + r\ — > J/ip + rj, and the corresponding 
cross section is only 6 nb at the same value of y/s. All other elastic processes listed in Table |l| are at the nb, pb, and 
even fb level. 

We now turn to the most important channels in our study: inelastic reactions. We show in Fig. || the isospin 
averaged total absorption cross section J/ip + tt — * D* + D plus D* + D. 
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FIG. 2. Isospin averaged total cross section for J ftp + n — > (D + D*) + (D* + D). 

Next in importance is the reaction J ftp + p — > D*D*, this is shown in Fig. ^. 
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FIG. 3. Cross section for J/ip + p -» D* + D*. 

Note that exothermic reactions are also possible, but they typically settle down to a low cross section value. A 
representative example is J/tp + p — > D + D, and this is shown in Fig. ^. 
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FIG. 4. Cross section for J/ip + p -> D + D. 



B. Phenomenological Model 

To compare with other calculations which have used effective Lagrangian methods but have constrained the model 
differently, we include this subsection. If we retreat somewhat from the symmetry and allow the coupling constants to 
be separately adjusted to empirical constraints or models we arrive at different predictions. For instance, if we choose 
the D*Dtt coupling constant to give a width consistent with a relativistic potential model prediction of 46 keV for 
D* ± Q, we find gD'D^ = 4.4 (whereas, the chiral prediction from the previous subsection used a value 3.02). Vector 
dominance arguments have been further used to fix couplings like gj/^p dd an d <?jv d* d* to be 7.7 Again, the 

chiral prediction used above is 4.93. We stress that the chiral model 
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FIG. 5. Cross sections for pion- and rho meson-induced dissociation of J/ip. The specific channels are the same as Figs, 
and []. 
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calculations are not different from the previous effective Lagrangian methods, practically speaking, the differences 
are merely methods of calibration |2lJ] . Note however that the calibration method we associated with the "phenomeno- 
logical model" will lead to K* phenomenology that is off by a factor of two, as seen in Table Gj, 

In Fig. |B| we show pion and rho dissociation of Jjtp in the phenomenological approach. The results are to be 
compared with the chiral model predictions in Figs, g and ^. So in some sense, the different results could be viewed 
as representative of uncertainties in the present hadronic approaches to J/tp dynamics. In this work the preference 
will go to the so-called chiral approach as it generates the hadronic phenomenology contained in Table | which does 
not contradict experimental measurements. 



III. ANOMALOUS PROCESSES 



We include next a section which reports on a couple of processes in the anomalous sector which turn out to give 
significant cross sections. Based on the observation that J/tp — > ij c +7 is 1.3% of 87 ke V, we use vector dominance to 
estimate the coupling of J/tp to rj c and uj. The Lagrangian we use is again listed in Eq. (2.13). We find a value g,j/^ c uj 
— 9.5 GeV -1 . Using similar reasoning and calculations for ui — > 7r°7, we extract g wp7r = 11-6 GeV -1 . Equipped with 
these couplings, we can compute the cross section for J/tp + ir^r] c + p through uj exchange. We present it in 
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FIG. 6. Cross sections for J/ip involving anomalous couplings and the r\ c meson. The specific channels are discussed in the 
text. 



Fig. and remark that it is quite large. Since this calculaton has been done, we have found that Shuryak and 
Teaney had considered this process previously using a model different from what is done here p^] . 

Similarly, the cross section for J/ijj + tt — > rj c + b\ is estimated. The value for the coupling gb l7 rw is deduced from 
the measured decay width of the b\ p6| . The cross section is displayed in Fig. [| Its value is also large. The elastic 
channel J/ip + co — > J/tp + uj proceeding through through rj c exchange can be considered. It too is included in Fig. @, 
and one can also conclude that it is relatively important. We have also calculated several other new reactions which 
we report on in the next section. 

With these cross sections, and the ones associated with all the other processes we have considered, the J/tp spectral 
function in a finite temperature gas of mesons can now be calculated. However before moving on to that topic, the 
important issue of form factors needs to be addressed. 



IV. HADRONIC FORM FACTORS 



The field theory in this work is formulated in a hadronic language and does not deal with fundamental fields but 
with degrees of freedom that are composite in terms of quark content. It is clear that the exchange of heavy mesons 
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(as the open charm D mesons, for example) leads to an interaction too short-ranged for the interacting particles to 
be left unmodified |25| l . Meson-exchanges are perhaps parameterizations of other phenomena which should be more 
evident at the parton level. This fact reveals itself in the appearance of hadronic form factors at the interaction 
vertices of Feynman diagrams. Those can be the source of additional uncertainty in the model. Note that form factor 
considerations are not restricted to meson exchange models like the one discussed here. For example, they occur also 
in the flux tube breaking model |24| and the 3 Po model |P5| . The choice of form factors is guided here by physical 
arguments, and they are introduced in a way that respects gauge invariance in the electromagnetic sector, and Lorcntz 
symmetry. 

A generic form is first chosen to construct t— and it— channel hadronic form factors. A candidate that lends itself 
to practical calculations is the monopole: 

A 2 

w A 2 + \ t-ml\ ' 

where m a is the mass of the exchanged meson. An advantage of this functional form is that the form factor is 
normalized to 1 for on-shell particles. Also, even if the kinematics venture into regions of time-like momentum 
transfer, this choice of form factors remains unitary. Each vertex therefore receives a contribution F(t) or F(u), 
depending on the appropriate kinematics. The vertices for the tadpole diagrams are determined by replacing its 
metric tensor structure by a general tensorial expansion constructed from the metric tensor and the available four- 
vectors. The coefficients of this expansion are then chosen such that the total amplitude is gauge invariant in the 
electromagnetic sector. The form factors in this work therefore build Feynman amplitudes that are both Lorentz- and 
gauge- invariant . 

Our values for A stem from elements of hadronic phenomenology which we describe now. We first fix the coupling 
constant gj/^ pT: by reproducing the measured width Tj/^^p^. The appropriate vector-vector-pseudoscalar (VV<p) 



Lagrangian has been shown in Eq. (2.13). For the determination of g the form factor (Eq. 4.1) plays no role, by 



construction. A is then determined by pushing one of the particles off-shell. Consider the measured pq ] radiative 
decay width T j/^^^o . Using a Vector Meson Dominance (VMD) argument, one may assume that the photon 
couples to the p of the above strong interaction vertex. Then, with g being determined, A can be obtained from a fit 
to the radiative decay width. A word of caution is necessary here: the J/tp also has a G-parity violating decay like 
J/?/ 1 — > W7r, so that presumably the photon could also originate from the to through VMD. However when compared 
with J/tp — > p°ir°, the decay into ton is suppressed by an order of magnitude so that we can safely ignore it here. 



In order to reproduce the J/tp radiative decay width one needs the parameter in Eq. 4.1 to be A = 1.25 GeV. This 
number is satisfying as it does represent a scale that is typically associated with soft hadronic interaction as are 
commonplace in, for example, the Bonn potential p7| . 

Another method to pin down hadronic form factors consists of considering J '/ip and open charm photoproduction 
data and to use their relation with J/?/;-nucleon total elastic and inelastic cross sections p8| . Using this argument, a 
J/tp + N inelastic cross section can be extracted from the data, and its value is «0.1 mb at i/s = 6 GeV [|29fl. Below 
this energy, some doubts have been expressed on the reliability of the cross section extraction through VMD |29| . We 
estimate the largest contribution to the J/ip-N inelastic cross section to be J/tpN — > A C D. Using the form factor 
described above and requiring that the inelastic cross section be 0.1 mb at y/s — 6 GeV sets A = 3 GeV. Past this 
energy value the J/ip-N meson-exchange cross section drops (unlike what is shown in Fig. (5) of Ref. pgfl), so that 
the upper bound set by photoproduction data is not exceeded. With this in mind A = 2 GeV is set as a conservative 
upper bound for the remainder of this work, thereby allowing for the contribution of other channels to the J/ip — N 
cross section. 

Summarizing, a range in A was estimated for the hadronic form factor introduced in the meson exchange model. 
Guided by hadronic phenomenology, we set 1.25 GeV < A < 2 GeV. The final cross sections are quite sensitive to 
the choice of the cutoff parameter A. This sensitivity is shown in Fig. u\ using the total inclusive absorption cross 
section for J/tp on tt. This cross section is the sum of the ones in Figs, g and Even within the window that has 
been determined for A, the cross section remains uncertain within an order of magnitude. The larger value of the 
form factor parameter has the total absorption cross section flattening out at around 4 mb. This value is dominated 
by the cross section in the anomalous sector with an r\ c in the final state, and is thus quite insensitive to the choice 
of the chiral model or the phenomenological model to calculate J/tp + tt — > D* + D + h.c. The apparent kink in 
the low energy region of Fig. |^ is related to the different thresholds for the reactions in Fig. || and |[ Also note that 
Shuryak and Teaney estimate that the cross section for J/tp + u — > r\ c + p is 1.2 mb, using non-relativistic quark 
model arguments |22| . The value obtained in the approach followed in this work is w 1 mb in the middle of the form 
factor range defined above. Those two different methods thus give numerical results that are not inconsistent with 
each other. It is important in this nonperturbative sector to cross-check model calculations. In this respect, the form 
factor corrected cross section for the reaction of Fig. || has a mean value of 0.1 mb, only slightly lower than that 
obtained in a quark-interchange model [p0[. 
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FIG. 7. The value of the total inclusive J/tp + n inelastic cross section, and its sensitivity to our choices of the form factor 
parameter A. 



The elastic cross section of the J/ifi with the w which was shown in Fig. g is also drastically affected by form factor 
considerations. This is displayed in Fig. ||. Note that the suppression due to the form factor is different in the elastic 
and inelastic cases, owing to different kinematics. 




VT (GeV) 

FIG. 8. The value of the J ftp + co elastic cross section, and its sensitivity to our choices of the form factor parameter A. 
The solid curve is obtained without form factors, the dashed curves mark the limits of our range in form factors. 
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FIG. 9. The value of the total inclusive J/ip + p inelastic cross section, and its sensitivity to our choices of the form factor 
parameter A. 



Finally we also show the total inclusive cross section for rho-induced absorption of the J /if) in Fig. It is worth 
noting that the upper choice of our form factor parameter has a large energy limit of « 1 mb. Using the methods 
just outlined, we have estimated also J/ip + r) — *■ r] c + (p (0.06 mb), J/ip + K — ► r\ c + K\ (0.18 mb). The numbers 
in parentheses refer to cross section values that are approximately in the centre of the form factor window. To our 
knowledge, those channels also have not been discussed previously. 



V. THE J/ip SPECTRAL FUNCTION 

For hadrons immersed in a strongly interacting medium, one finds that two-loop effects dominate over one-loop 
effects in the calculations of imaginary parts of particle propagators ||l|] . This owes largely to the size of the coupling 
in the confined sector. We have verified that this is the case here also, by doing explicit calculations. Thus, we neglect 
one-loop effects. We will also neglect effects on the real part of the J/ip propagator. We thus will assume that the J/ip 
will suffer negligible mass-shifts. We partly base this reasoning on the large mass of the charmonium vector meson. 
Also, recent calculations of J/ip properties in nuclear matter do yield mass shifts that are small Ej^ |. 

Our first task is then to calculate the broadening due to collisions of the J / ip with particles that make up the heat 
bath. The width induced by a reaction of the type J/ip 2 — > 34, where 2, 3, and 4 are arbitrary species is | gjl[j3| . 

1 



r(w,fl = — / <mn 2 (E 2 )(l + n 3 (E 3 ))(l + r H (E 4 ))\M(J/ib2^M)\ 2 , (5.1) 



where lu = Jp 2 + m 2 ,,, p being the three- vector of the J/ip. Note that 3 or 4 can be a J/ip. In Eq. (5.1), 



V/V" 

d£l = dp 2 dp 3 dp i {2TT) 4 5(p + p 2 — Pz — Pi), (5-2) 

and 

m = ■ (5 - 3) 

The reactions we consider will involve only mesons. In principle, a J/ip can also be produced by an inverse reaction 

involving particles from the thermal background J33j. Here, phase space considerations make the inverse channel 
negligible. 
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We write the spectral function of the J/i[> as 

Aj/^(uj,p) = -2ImD J/V) (w,p) , 



(5.4) 



where Dj/^ is the scalar part of the J/ip propagator. Neglecting the difference between longitudinal and transverse 
polarizations of the J/ip in the finite temperature medium [j34f , one has 



Dj/^{u),p) 



P 2 ~ m J/^ ~ F ( w iP) 
where p^ — (u),p), and F is the scalar imaginary self-energy. Then, using 



Tj N , = -1/mj/^ Im F(p 2 = m 
one can write in the on-shell approximation 



J/4, 



A 



(p 2 



m 2 , ,) 2 + m 2 , l T 2 



(5.5) 



(5.6) 



(5.7) 



where T 



J/ip 



contains all the contributions we have discussed so far: the vacuum width and the contributions from 



clastic and inelastic collisions. 

One can first shown the effects of purely elastic processes on the J/ip spectral function. This appears in Fig. [l(]. 
Note that A = 2 Gev throughout this section, keeping in mind any investigation of a more quantitative nature will 
need to reflect the possible ambiguities in this choice that were discussed earlier. 
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FIG. 10. Spectral function in vacuum and at finite temperature allowing only elastic scattering. 



Without inclusion of the elastic scattering with w's through rj c exchange, the spectral function deviates so little 
from the vacuum that all curves lie on top of one another. But as one can see in Fig.|lO|, elastic scattering in the 
medium now has a non negligible effect and distorts the spectral function. Since the J/ip + lu cross section with form 
factors is quite flat in y/s, the momentum dependence here is not important. 

Let us now consider the J/ip charmonium state traveling in a finite temperature gas first consisting only of 7r's, K's 
and p's. The spectral function at two temperatures, 150 and 200 MeV, is shown in Figs, [ll] and [l2]. One notices a 
substantial broadening of the spectral distribution, along with a suppression of the peak. This considerable effect is 
even more pronounced at the higher temperature. If we include all the inelastic processes we have considered in this 
work (summarized in TableO, it turns out the quantitative differences between results involving all those and the 
ones shown in Fig. 11 and [l2; are small and can be neglected. Clearly, the 7r's, K's, and p's play the leading role. 
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FIG. 12. Spectral function in vacuum and at 200 MeV temperature allowing inelastic interactions with it's, K's and p's. 



VI. CONCLUSION 



In this work, cross sections for the interactions of the J/ip with light mesons were evaluated. It has been found that 
those cross sections are all quite different, and also not constant with respect to the energy of the colliding particles. 
The form factors that are germane to meson exchange models such as the one discussed here have been constrained 
by hadronic phenomenology, Lorentz invariance and electromagnetic gauge invariance. The numerical effect of those 
form factors are large and therefore a careful treatment is mandatory. The importance of the J/ip + n — ► X hadronic 
absorption channel has been highlighted, especially those that invole the r] c . Similarly, the elastic cross section 
involving an external u interacting through an exchanged rj c is found to be appreciable. Those findings should have 
some effect on the hadronic phenomenology that is an ingredient to the theoretical modeling of high energy heavy 
ion collisions. Specifically, the J/ip + 7r inelastic channel can be interpreted as the leading contribution in the J/ip 
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absorption by "comovers". In this light, the average value of our comover cross section just about reaches 2 mb. If 
one consider the J/ip + p channel, one could even add an additional mb. A recent study uses cr comovors = 1 mb |J. It is 
important for values obtained phenomenologically and for those based on more microscopic approaches to eventually 
meet. The studies performed here however do point to the richness of the many-body problem. Several reactions 
channels have been considered and even more work needs to be done in order to complete the survey of what turns 
out to be a vast hadronic landscape. Our exploration continues. 

We have evaluated the spectral function for a J/ip state traveling in a finite temperature gas of mesons. We 
have found that the spectral function gets considerably modified, owing mainly to inelastic interactions with the 
constituents of the hot meson gas. 

For the moment, we have refrained from attempting a detailed comparison with heavy ion data. It is also clear 
that such an application will require great care. For example it is claimed here that the absorption of the J/tp on 
pions is important, especially when the r\ c is part of the final state. However, inverse reactions can produce a J/ip'. 
Vc + P - * J/ip + ir (0.7 mb), and r] c + Tr — > J/ip + p (2 mb). The numbers in parentheses refer to cross section values that 
are approximately in the centre of the form factor window. One therefore needs to advocate a careful simulation of the 
nuclear collision, with the inclusion of the relevant important reaction channels and of their detailed balance partners. 
This work is extensive, has begun, and will be reported on elsewhere. It is first necessary to place the interactions 
of the J/tp in a proper many-body setting, at the appropriate level of sophistication our current understanding of 
hadronic physics requires. This also implies pointing out the caveats as well as the successes. 
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APPENDIX: GAUGE INVARIANCE 



The amplitudes discussed in this work that couple to vector mesons that have the quantum number of the photon 
need to obey gauge invariance in the electromagnetic sector, or more specifically current conservation. This statement 
is an immediate consequence of Vector Meson Dominance. We focus for the moment on the reaction J/V'(pi) + i"(P2) — > 
D*(j>a) + D(pi). The invariant amplitude emerging from Eq. 2.11 is M = Mx + M2 + M3 where 



Mi 
M 2 

Ma 



|g 2 e ^( Pl )(2 P 4 - Pi), {pi _ p ] )2 _ ml (2P2 - PsU»(pa) 



-^9 2 ^{pi) [(2p3 -pi)»gva - (px+pa) a g^u + (2pi - p a ) v g^] 

[-9 a0 + (pi-p 3 ) a {pi-p 3 ) p /ml] / 
x (px-^-mi) (Pz+P^* (P») 

f5 2 ^(pi)[-^]e"(p a ). 



(Al) 



Gauge invariance requires that M (e^(pi) — > Pi) must be identically zero. In particular, it must vanish for arbitrary 
choices of pseudoscalar and vector masses. 

Replacing e M (pi) — > and doing the contraction gives 



Mi 
M 2 



g 2 [-(2p 2 -p 3 ) v ] e> 3 ) 

\g 2 \(P2 +Pi)y + 



(Pi -P3) 2 ~ ml 



(2^3 -Px) -Pxipx - Pa) ■ O2 +Pi)(Px ~Pa)v 



+ (PX +P3) ■ (P2 + Pi)(Px)v ~ 



(pi +P3) • (pi ~Pa)(px - Pa) ■ (P2 +P4)(Pi)>. 



M? 



(2pi -P3)w(Pi) • (P2 +Pa) + 

V [-(Pi),]e> 3 ). 



(2pi -Pa)v(px) ■ (Pi -P3)(pi -P3) • (P2 +Pi) 



C"(P8) 



(A2) 



We note that A^i plus the first term in M2, plus M3 vanishes due to energy- momentum conservation. The 
remaining pieces from M2 are 



1 



oc 



-{{2pa -px)-x (Pi ~Pa) • (pa +P4)(pi -ps)^ + mf(pi +P3) ■ (P2+Px)(px)i 



m 3 [(Pi -P3) 2 - ml] 
-{Px +P3) ■ (Pi -P3XP1 -P3) • (P2 + Px)(Px)» - rnl(2pi - pa)uPx • (P2 +P4) 

+ (2pi -P3)vP\ ■ (Pi ~P3)(Pl ~P3) ' (P2 +P4)}e' y (P3) 



(A3) 



Terms proportional to (pa)u vanish when contracted with t v (pa) due to transversality. Thus, the surviving terms 
can be written as 



x 



(Px\ 



:{{2pa ~Px) ■ (Pl)(P4 -P2) • (P2 +P4) +P3) • (P2 +P4) 



m 3 [(PX - P3) 2 - ™%] 

+^l(-2Pl) ' (P2 +P4) - (Pi +P3) • (Pi -P3)(P4 -P2) ' (P2 +P4) + 2(pi) • (pi -p 3 )(P4 -P2) ' (P2 +Pi)}( U (A4) 

Finally, we can simplify to 

oc (px)u{(m\ - ml) [ml] + m 2 [(pa - pi) ■ (P2 +PA)]}e v (pa) 



= (px)v{ml(ml - ml) +ml(p 2 - Pa) ■ (P2 +Pa)}^(p3) 
= (pi) v {m\{m\ - ml + ml - m 2 } 
= 0. 



(A5) 



Indeed, current is conserved in the general case. The other channels can similarly be shown to conserve current. 
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TABLE I. Model prediction for widths. The phenomenological approach refers to that of Refs. |ll|Jl 



particle 


chiral model 


pheno. model 


experiment 


K*(892) u 


44.5 MeV 


97.0 MeV 


50.5 ± 0.6 MeV 


K*(892)± 


44.5 MeV 


97.0 MeV 


49.8 ± 0.8 MeV 


D*(2007) u 


10.1 keV 


22.0 KeV 


< 2.1 MeV, 90% CL 


D*(2010) ± 


21.1 keV 


46.0 KeV 


< 131 keV, 90% CL 



TABLE II. We list here the hadronic reactions involving J/tp that were considered in this work. It is implied that the figures 
also include the Hermitian conjugate inelastic final state when it is different from the one listed below. 



iijlastic cnannels 


initial state 


Inelastic channels 


final state 


J/lp + TV 


J/lp + 7T 




D* + D 


J/tp + V 


J/4> + TT 




ric + p 


J/^ + K 


J/lp + TT 




Vc + h 


J/i> + P 


J/4> + v 




D* +D 


J/lp + LU 


J/ip + r) 




rjc + <P 


J/if> + 4> 


J/ip + K 




D s +D* 




J/4> + P 




D + D 




J/iP + P 




D* + D* 




J/tP + P 




Vc + 1" 




J/ip + u 




D + D 




J/ip + u 




D* + D* 




J/ip + 4> 




D + D 




J/4> + (p 




D* +D* 




J/iP + K 




rjc + K 1 




J/%P + K* 




D s +D 




J/ip + K* 




D* + D* 



1G 



